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Abstract. We review the expected science performance of the New 
Gravitational- Wave Observatory (NGO, a.k.a., elJSA), a mission under study 
by the European Space Agency for launch in the early 2020s, eUSA will sur- 
vey the low-frequency gravitational- wave sky (from 0.1 mHz to 1 Ha), detecting 
and characterizing a broad variety of systems and events throughout the Uni- 
verse, including the coalescences of massive black holes brought toget her by galaxy 
mergers; the inspirals of stellar-mass black holes and compact, stars into central 
galactic black holes; several millions of uEtracompact binaries, both detached and 
mass transferring, in the Galaxy; and possibly unforeseen sources such as the 
relic gravitational-wave radiation from the early Universe, eLISA’s high signal- 
Lonoke measurements will provide new insight into the structure and history of 
tht* Universe, and they will lest general relativity in its strong- field dynamical 
regime. 

PA OS numbers; 04,80. Nu, 95,55, Ym 

1* Introduction 

Over the lost two decades, as many &$ 2,500 articles on space-based gravitational- 
wave (GW) detection Included mentions of LISA (the Laser Interferometer Space 
Antenna) [1, 2, 3j, the space-based GW interferometer planned and developed together 
by NASA anti ESA. This collaboration between the two agencies ended in early 2011 
for programmatic and budgetary reasons. In fact, LISA, as brought forth by the 
entirety of those papers, was more than a space project: it was the concept (and 
the cherished dream) of a space-based GW observatory that would explore the low- 
frequency GW sky, in a frequency band (l0“ 4 1 Hz) populated by millions of sources 
in the Galaxy and beyond: compact Galactic binaries; coalescing massive black holes 
(MBHs) throughout the Universe; the captures of stellar remnants into MBHs; and 
possibly relic radiation from the early Universe, 

Ah along its evolution, the LISA design remained based on three architectural 
principles developed and refined since the 1970s: a triangular spacecraft formation 
with Mkm arms, in Earth-like orbit around the Sun: the continuous monitoring of 
inter-spacecraft distance oscillations by laser interferometry; drag-free control of the 
spacecraft around freely falling test masses, the reference endpoints for the distance 
measurements, achieved using micro-Newton thrusters. The current incarnation of 
this concept is eLISA (evolved LISA), a mission under consideration by ESA alone 
(under the official name of NGO, the New Gravitational-wave Observatory) for launch 
in 2022 within the Cosmic Vision program. 

The eLISA design would achieve a great part of Hie LISA science goals, as 
presented In [If, and endorsed by the 2010 U.S. astronomy and astrophysics decadal 
[4], This article surveys eLISA s science performance (sensitivity, event rates, and 
parameter estimation), as scoped out by these authors in the spring and summer of 
2011, and as discussed in full in Ref. [5j. This article is organized as follows: in Sec, 2 
we provide a very brief overview of eLISA and its GW sensitivity while later sections 
are organized by science topics. In Sec. 3, we discuss the astrophysics of compact 
stellar-mass binaries in the Galaxy: in Sec. 4, the origin and evolution of the massive 
BHs found at the center of galaxies, as studied through their coalescence GWs: in See. 
5, the dynamics and populations of galactic nuclei, as probed through the captures of 
stellar-mass objects into massive BHs: in Sec, fi. rhe fundamental theory of gravitation, 
including its behavior in the strong nonlinear regime, its possible deviations from 
general- relativistic predictions, and the nature of BHs; in Sec. 7, rhe (potentially new) 
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physics of r hi* early Universe, and the measurement of cosmological parameters with 
GW events- Lust, in Sec, 8 we draw our conclusions, and express a wish. 

2. The eLISA mission and sensitivity 

We refer the reader to )5] for a detailed description of the eLISA architecture. eLISA 
lias a clear LISA heritage, with a few substantial differences. The eLISA arms 
will be shorter (1 Mkm), simplifying the tracking of distant spacecraft, alleviating 
requirements on lasers and optics, and reducing the mass of the propellant needed 
to reach the final spacecraft orbits. The orbits themselves may be slowly drifting 
away from Earth, again saving propellant, and the nominal mission duration will be 
two years, extendable to five. As much existing hardware as possible, including the 
spacecraft bus, will be incorporated from the LISA Pathfinder mission, scheduled for 
launch by ESA in 2014. Tile three spacecraft will consist of one “mother" and two 
simpler “daughters," with interferometric measurements along only two arms, for cost 
and weight savings that make launch possible with smaller rockets than LISA. (Note 
that LISA was to be built with laser links along the three arms, but it was not a 
requirement that they would operate throughout the mission.) 

The eLISA power-spectral-de nsity requirement for the residual test-mass 
acceleration is S acc (f) = 2.13 x I0 _ ^{ 1 + 10 -4 Hz//) nr s -4 Hz - 1 . while the position- 
noise requirement breaks up into 5 S11 {/} = 5.25 x IQ -23 m 3 Hz -) for shot noise, 
and S omn (f) - 6.28 x 10“ 23 m 2 Hz” 1 for all other measurement noises. With 
these requirements, eLISA achieves the equivalent-strain noise plotted in Fig. 1. and 
ai ) p roxim a t v< I a nalyt teal ly by 

qt r\ _ 4 S HVC (f) T n{f) + S otm i(f) _ / . / / \^\ , 

&(/) "T £2 x \ f 0.41 c/2£/ / ’ (1) 

where L 1 Mkm. c is the speed of light, and S(f) has already been normalized to 
account for the sky-averaged eLISA response to GWs. This noise allows the detection 
of a strain of about 3.7 x 10~ 24 in a two-year measurement with an SNR of 1, The 
requirement on the useful measurement band is 10“ 4 Hz to 1 Hz. with a goal of 3x 10“ s 
Hz to 1 Hz. 

3, Compact binaries in the Galaxy 

The most numerous sources in the low-frequency GW sky observed by eLISA will 
be short-period binaries of two compact objects such as white dwarfs (WDh) or 
neutron stars (NSs). These systems have* weak GW emission relative to the much 
heavier massive-BH binaries, but are numerous in the Galaxy and even in the Solar 
neighborhood. To date, ast ronomers have observed about 50 ultra-compact binaries 
with periods shorter than one hour, comprising both detached systems and interacting 
binaries where mass is being transferred from one star to the other. Wide-field and 
synoptic surveys such as SDSS and PTF (and in the future, PanSTARRS, EG APS, and 
LSST) will continue to enlarge this sample 7. 8]. Interacting ultra-compact binaries 
with NS accretors are found by all-sky X-ray monitors and in dedicated surveys [9], 

A large subset of known systems will be guaranteed verification sources for 
eLISA [lOj; their well-modeled GW signals will be detected within the first few 
weeks to months of operation, verifying instrument performance. The most promising 
verification binaries are the shortest-known- period interacting systems HM Cue (with 



Low-frequency gravitational- wave science with cLISA/NGO 


4 



Figure J. eLISA equivalent-strain noise, averogi^l over source sky location and 
polarization, as a function of frequency. The solid red curve was obtained wiib 
the LISACode 2.0 simulator [<ij, while the dashed blue curve is plotted from Eq. 
(1). For comparison, the dotted green curve shows the LISA sensitivity. 


a period of 5.4 m [11] }, V407 Vul (P 9.5 m), and ES Get and the recently discovered 
detached system SDSS J065I+28 (P = 12 in; [12]). 

eLlSA will individually detect and determine the periods of several thousand 
currently unknown compact binaries (in our estimate, 3,500-4,100 systems 
for a two-year observation; [5, 14]), wliile the combined signals of tens of millions 
unresolvable systems will form a stochastic GW foreground at frequencies below a 
few mHz ((15, 16]; see Fig. 2.) About ^ 500 dose or high-frequency (> 10 mHz) sources 
will be seen with large SNRs. allowing the determination of sky position to better than 
10 deg 2 , of frequency derivative to 10%. of inclination to 10 deg, and of distance to 
10%. This large sample will allow a detailed study of the Galactic population* which 
is poorly constrained by EM observations and theoretical predictions [17]. 

Detections will be dominated by double WD binaries with the shortest periods 
(5 L 0 minu tes ) , Their i nergers are candid ate ] ) rogenitors for i na i ly i i iterest ing sy st ems: 
type la [18] and peculiar supernovae [19. 20]: single subdwarf O and B stars, R Corona 
Borealis stars and maybe all massive WDs [21); and possibly the rapidly spinning NSs 
observed as ms radio pulsars and magnet ars [221 These binaries are short lived, very 
faint for telescopes* and scarce (few thousand in the whole Galaxy), so GWs will 
provide a unique window on their physics. eLISA will determine their merger rate, 
constrain their formation, ami illuminate the preceding phases of binary evolution, 
most notably the common-envelope phase. 

Common-envelope evolution is crucial to most binary systems that produce 
high-energy phenomena such as "-ray bursts and X-ray emission, but our 
understanding of its physics and outcome is limited (23, 24] and challenged by 
observations [25, 26], The standard scenario is as follows. Most shirs in the Universe 
arc in binaries, and roughly half of binaries are formed at close enough separations that 
the stars will interact as they evolve into giants or supergiants, Following runaway 
mass transfer, the companion of the giant can end up inside the outer layers (the 
envelope) of the giant; dynamical friction reduces the velocity of the companion. 
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Figure 2. Main figure: power spectral density o t the stochastic GW foreground 
from Galactic binaries, before (blue) and after (red) the subtraction of individually 
resolvable systems, which are plotted ns green and red/blue dots (for detached and 
mass-transferring systems), A few known verification binaries are shown as white 
dots* The sol id /dashed black curves trace instrument noise alone/with confusion 
noise. Spectra are shown fortheTDI observable "X 4 ’ (see, e.g., [ 13}) ; subtraction is 
simulated for a two-year observation and threshold SNR = 7: resolvable systems 
are placed a factor SNR 2 above the combined instrument and confusion noise. 
Inset: time series of the residual foreground, which carries information about the 
number and distribution of binaries in the Galaxy. 


shrinking the orbit and transferring angular momentum and energy into the envelope; 
the envelope eventually becomes unbound, leading to a very compact binary consisting 
of the core of the giant and the original companion [27]* 

eLISA will also test dynamical interactions in globular clusters, which produce 
an overabundance of ultra-compact X-ray binaries consisting of a NS accreting 
material from a WD companion. The eLISA angular resolution will he sufficient 
to distinguish WD binaries in clusters, verifying whether they are also plentiful. 

The eLISA measurements of individual short-period binaries will provide a wealth 
of information on the physics of tidal interactions and the stability of mass transfer* 
For detached systems with little or no interaction, t he evolution of the GW signal is 
dominated by gravitational radiation: 

h xM* /3 f' 2/ *D- 1 , f oc M 5f3 f n t 3 , /==/. ( 2 ) 

where h is the GW strain, / the GW frequency, + m 2 } 1 ' ' is 

the chirp mass with m 1 . m? the individual masses, and D is the distance. Thus, 
measuring ft, /, and / {which will be possible in 25% of systems) provide M and D\ 
measuring also / (which may be possible for a few high-SNR systems) tests secular 
effects from tidal and mass- transfer interactions. Short-term variations are not likely 
to prevent detection }2S], and the precision of / and / determination increases with 
the duration of the mission. 
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Tidal interactions are possible when at least one binary component does not 
corotate with the orbital motion, or when the orbit is eccentric. Their strength is 
unknown [29], and has important consequences on the tidal heating (and possibly 
optical observability) of WD binaries, as well as the stability of mass transfer. This 
process begins after gravitational radiation shrinks detached binaries to sufficiently 
close orbits {with P ^ a few minutes) that one of the stars fills its Roche lobe and 
its material can leak to the companion. Mass transfer can be self- limiting, stable, or 
unstable, depending on the resulting evolution of the orbit and of the donor radius. 
Unstable transfer leads to mergers; stable systems (the interacting WD binaries known 
as AM CVn systems, as well as ultra-compact X-ray binaries) will be observed - and 
counted by eLISA in the early stages of mass transfer [30]. Efficient tidal coupling 
can return angular momentum from the accreted material to the orbit [29, 31, 32) , 
slowing the inspiral and increasing the fraction of WD binaries that survive the onset 
of mass transfer from 0.2% to 20% [33] . 

The unresolved foreground from Galactic binaries will provide a light 
additional noise component for the detection of loud broadband signals (see Fig. 2), 
but it also contains precious astrophysieul information. Its overall level measures 
the total number of binaries (mostly double WDs); its spectral shape characterizes 
their history and evolution: and its yearly modulation [34], together with the distance 
determinations from many individual systems, constrains the distribution of sources 
in tfie different Galactic components. Thus eLISA will probe dynamical effects in the 
Galactic center, which may increase the number of tight binaries [35]; it will measure 
the poorly known scale height of the disk; and it will sample the population of the 
halo [36, 15], which hosts two anomalous AM CVn systems and which may have a 
rather different compact-binary population than the rest of the Galaxy Furthermore, 
the eLISA measurements of orbital inclinations for individual binaries, compared with 
the overall angular momentum of the Galaxy, will provide hints on the formation of 
binaries from interstellar clouds. 

eLISA will also constrain the formation rate and numbers of NS binaries and 
ultra-compact stellar- mass BH binaries, throughout the Galaxy and without EM 
selection effects. These numbers are highly uncertain, but as many as several tens of 
systems may be detectable by eLISA [33, 37], complementing the ground-based GW 
observations of these same systems in other galaxies (and at much shorter periods). 

More generally the astrophysieal populations and parameters probed by eLISA 
will be different from, and complementary to, what can be deduced from EM 
observations. For instance. eLISA will be sensitive to binaries at the Galactic center 
and throughout the Galaxy while GAIA [38] will be limited to the Solar neighborhood: 
GWs encode distances and orbital inclinations, while EM emission is sensitive to 
surface processes. Dedicated observing programs and public data releases will allow 
simultaneous and follow-up EM observations of binaries identified by eLISA. 

4. Massive black-hole binaries 

(See [5: for a much deeper review.) According to the accretion paradigm [39, 40. 41], 
snpermasslve BHs of 10 s 10 4 * * * * 9 4/. power quasars active galactic nuclei so luminous 

that they often outshine their galaxy host, which are detected over the entire cosmic 

time accessible to our telescopes. Quiet superm&ssive BHs are ubiquitous in our low- 

reds hift Universe, where they are observed to have masses closely correlated with key 

properties of their galactic host (see [42], and refs, therein) leading to the notion that 
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galaxies and their nuclear MBHs form and evolve in symbiosis (see, e.g., [43, 44, 45])* 

In the currently favored cosmological paradigm, regions of higher-density cold 
dark matter in the early Universe form self-gravitating halos, which grow through 
mergers with other halos and accretion of surrounding matter; baryons mid MBHs are 
thought to follow a similar bottom-up hierarchical clustering process [46, 47, 48* 49, 
50], MBHs may be born as small seeds (10 2 10 : * M ) from the core collapse of the 
first generation of "Pop III” stars formed horn gas clouds In light halos at £ — 15-20 
[51, 50 j; or as large seeds (10 : * 10 5 Mq) from the collapse of very massive quasi-stars 
formed in much heavier halos at z * 10-15 [52, 53]; or by runaway collisions in 
star clusters [54}; or again by direct gas collapse in mergers [55] (See [56, 57] and 
refs, therein). The seeds then evolve over cosmic time through intermittent, copious 
accretion and through mergers with other MBHs after the merger of their galaxies. 

The cosmic X-ray background from active MBHs at z < 3 suggests that 
radiatively efficient accretion played a large part in building up MBH mass [58, 59, 60], 
so information about the initial mass distribution is not readily accessible in the local 
Universe. By contrast, eLISA will measure the masses of the original seeds From 
their merger events. Furthermore, ii is unknown [6.1] whether accretion proceeds 
coherently from a geometrically thin, corotating disk [62] (which can spin MBHs up 
to the J/M 2 0.93-0,99 limit imposed by basic physics [63, 64j) or chaotically from 

randomly oriented episodes [65] (which typically result in smaller spins). eLISA ‘s 
accurate measurements of MBH spins will provide evidence for either mechanism [66]. 

After a galactic merger, the central MBHs spiral inward, together with their 
bulge* or disc, under the action of dynamical friction, and pair as a pc-snde Kepleriau 
binary [67, 68, 69, 70, 71]; MBH binaries are then thought to hanlen into gravitational- 
radiation- dominated systems by ejecting nearby stars (assuming a sufficient supply) 
[72, 73, 74] or by gas torques and flows in gas-rich environments [75, 76, 77]; the final 
binary coalescence is the most luminous event in the Universt 1 (albeit in GWs). BH 
mergers have been explored only recently by numerical relativity [78] . showing how the 
mass and spin of the final BH remnant arise from those of the binary components, and 
predicting remarkable physical phenomena such as large remnant recoils for peculiar 
spin configurations [79]. The predicted coalescence rate in the eLISA frequency band 
ranges from a handful up to few hundred events per year, depending on theoretical 
assumptions ( [80. 81, 82, 83, 84, 85, 86, 87]), 

eLISA will be sensitive to GW signals from all three phases of MBH coalescence 
(inspirah merger* and ring-down [88]). To assess the eLISA science performance in 
this area, after experimenting with different waveform families, wc modeled these 
signals with the "PhenomC” phenomenological waveforms [89], which stitch together 
post-Newtonian (PN) inspiral waves [90] with frequency- domain fits to numerically 
modeled late-inspiral and ringdown waves. (PhenomC waveforms were derived for 
binaries with aligned spins; we adapt them to the non-aligned case by projecting the 
orbital angular momentum and individual spins onto the angular momentum of the 
distorted BH after merger.) 

The first metric of performance is the detection SNR, angle-averaged over sky 
position and source orientation, which is plotted in Fig. 3 as a function of total 
rest mass and cosmological redshift (left panel) and as a function of total rest mass 
and mass ratio for binaries at z 4 (right panel). eLISA covers almost all the 
mass redshift parameter space of MBH astrophysics: any equal-mass biliary with 
M io t I0 +i 10 7 M. (the crucial “middleweight** range inaccessible to EM observations 
beyond the local Universe) can be detected (with SNR > 10) out to the highest 
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Figure 3. Left: constant -IfveJ contours of sky- and polarization-averaged SNR 
for equal -mass non-spinning binaries as a function of total rest mass JUtot and 
cosmological redahlft z. The 8NR includes in spiral, merger and ringdown. Right: 
SNR contours as a fu nr lion of Mtot and mass ratio g — mi/wtj. 




Figure i. Left: distribution of expected SNR for MBH mergers as a Func tion 
of z, computed from the SE/LE/SC/LC metacatalog (see main text). Right: 
likelihood for the mixing fraction for an individual realization of mixed model 
*FSE + (1 - ,F)LE wiih J 7 ~ 0-45 (see main text). 


redshifts, while equal-mass binaries with M\ ni > 10 5 M * are seen in detail as strong 
signals (SNR > 100) out to z 5. Binaries with M Uyi > 10 5 M® and mass ratios < 10 
are seen wit li SNR > 20 out to z = 4. 

To evaluate expected SNRs in tin 1 context of realistic MBH populations, we 
consider four fiducial scenarios (SE, LE, SC. LG) where MBHs originally form from 
Small \mM: ) or Large seeds 11) J A/, )♦ and where they subsequently grow by 
Extended or Chaotic accretion. (See [91] for details; here we enhance that analysis 
by including random spin orbit misalignments up to 20 deg in E models |92’). For 
each scenario we generate multiple catalogs of merger events, and join them in equal 
proportions into a single metacatalog. Figure 4 shows the resulting distribution of 
SNR with eLISA will detect sources wit El SXR > 10 out to z < 10* a limit imposed 
by masses of the expected binary population as a function of 

For the same metacatalog, Fig, 5 shows the expected accuracy of parameter 
determination, estimated using a Fisher-matrix approach based on PN inspiral 
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Figure 5 , Parameter-estimation accuracy (relative frequency of fractional or 
absolute errors over SE/LE/SC/LC metacatalog) For primary and secondary 
rvdshijled MBH masses and dimensionless spins (mi and m2* ni fm\ and 03/7712 * 
respectively )i luminosity distance D t and sky position AH* 

waveforms with spin-induced precession, augmented with PhenomC merger-ringdown 
waveforms* eLISA can determine the tndshifted component masses (m T ed shift — 
(1 -fz) m r <*t ) to 0.1-1 %; the primary-MBH spins to 0*01-0.1; and the secondary-MBH 
spins to 0.1 in a fraction of systems. (Compare with EM MBH-mass uncertainties ~ 
15-200%, except for the Milky Way MBH, and with very large MBH-spin uncertaint ies 
from Ka iron line fits [93].} The errors hi Dt have a wider spread, from a few percent 
to virtual non-determination, while sky posit ion Q is typically determined to 10 1000 
deg 2 , Compared to previous published estimates for LLSA, the accuracy in determining 
both Dt and il is reduced for eLISA by having interferometric measurements only 
along two arms (although three arms were always a goal, not a requirement, for USA). 

The next order of analysis is to combine multiple MBH-coalescence observations, 
resulting in a catalog of binary/remnant parameters, into a single inference about 
the mechanisms of MBH formation and evolution throughout cosmic history. 
This problem was analyzed extensively by Sesana and colleagues [94] in the context 
of LISA. We repeated their analysis for eLISA, by generating 1,000 catalogs of 
detected mergers (over two years) for each of the four SE/LE/SC/LC scenarios, 
and comparing the relative likelihood p(Avs, B) = p(A|C)/lp(A|C)H-p{B|C)] for each 
pair of scenarios ( A, B), for C — A or B, We considered only detections with SNR > 8. 
and used spinless, restricted PN waveforms. Table l shows our results for a relative 
likelihood threshold 0.95: for instance, the first row on the left shows that if SE is 
true, it cautd be discriminated from LE and LC in 99% of realizations, but from SC 
only in 48% of realizations; the last row on the left shows that LC could not hr ruled 
out in 2% of realizations when SE or SC are true, but in 22%i of realizations when LE 
is true. This degeneracy between accretion mechanisms is an artifact of the spin- less 
assumption: including information about the spin of the final merged MBH (which 
can be measured in 30% of detections) provides essentially perfect discrimination. 

Last, because no theoretical model will exactly capture the "true 1 ' formation and 
evolution history of MBHs. we investigated eLISA s ability of measuring the mixing 
fraction 0 < F < I in a mixture model FA + (I — F)B that produces coalescence 
events with probability F from scenario A, and 1 — T from B . For instance, for the 
ease .FSE -h 1 1 — J r )LE with F — 0 . 45 , F can be measured with an uncertainty of 
0.1 (see right panel of Fig. 4). Although highly idealized, this example shows the 1 
potential of eLISA ? s observations to constrain MBH astrophysics along their entire 
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Tnhlf; I, Model discrimination with fiLISA MBH-binary observations. The 
upper-righl half <jf each table shows t he fraction of realizations in which the row 
model would be chosen over the column model with a likelihood threshold > 0-95, 
when the row model is true. The lower deft half of each table shows the fraction 
of realizations in which the coup model cannot be ruled out against the column 
model when the column model is true. In the left table we consider only the 
measured masses and retishift for observed events; in the right table we Include 
also the observed distribution of remnant spins. 
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cosmic history, in mass and red&hift ranges inaccessible to EM astronomy. 

In closing this section, we note that eLISA may also detect coalescences of BHs 
with masses of 10 2 TO 4 A/ 0 (intermediate-ma-ss BHs, or IMBHs). These events do not 
result from hierarchical galaxy mergers, but they occur locally under the extreme 
conditions of star clusters. IMBHs may form in young clusters by way of mass 
segregation followed by runaway mergers [95, 96, 97, 98, 99]; IMBH binaries may form 
in situ [100], or after the' collision of two clusters [101, 102]. Although the evidence 
for IMBHs is tentative [103, 104], eLISA may observe as many as a few coalescences 
per year [10 1] out to a few Gpc [89]; it may also detect stellar-mass BHs plunging into 
IMBHs in the local Universe [105], 

5. Extreme-mass- ratio inspirals and the astrophysics of dense stellar 
systems 

There is of course one galactic nucleus, our own, that can be studied and imaged 
in great detail [106, 107, 108, 109, 110, 111], The central few parsecs of the Milky 
Way host a flense, luminous star cluster centered around the extremely compact radio 
source SgrAL The increase in stellar velocities toward SgrA* indicates the presence of 
a (4 ±0.4) x ID 0 Mv, central dark mass [ill]* while the highly eccentric, low- peri apsis 
orbit of young star S2 requires a central- mass density > 10 13 A/f. pe -:i [112]; a density 
> 10 13 M , pc * is also inferred from the compactness of the radio source [113]. These 
limits provide compelling evidence that the dark point-mass at SgrA* is an MBH 
[112, 114, 115], 

Unfortunately, the nearest external galaxy is 11)1) times farther from Earth than 
SgrA", and the nearest quasar is 100,090 times farther, so Imaging other galactic 
centers is prohibitive. It will however become possible with eLISA. This is because 
MBHs are surrounded by a variety of stellar populations, including compact stellar 
remnants (stellar BHs, NSs, and WDs) that can reach very relativistic orbits around 
i he MBH without being tidaily disrupted [ 110]. The compact stars may plunge directly 
into the event horizon of the MBH: or they may spiral in gradually while emit ting G Ws. 
These latter systems, known as extreme-mass ratio inspirals (EMRIs). will produce 
signals detectable by eLISA for MBH masses of 10 4 10' A/. . Stellar-mass BHs should 
be concentrated in cusps near MBHs 117. 11 8. 98. 119. 120 and generate stronger 
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GWs thanks to their relatively larger mass, so they will provide most detect! oils. 

EMRIk are produced when compact stars in the inner 0.01 pc of galactic nuclei 
are repeatedly scattered by other stars into highly eccentric orbits where gravitational 
radiation takes over their evolution [116]; resonant relaxation caused by long-term 
torques between orbits increases the rate of orbit diffusion [121, 122] t although 
relativistic precession can hinder this mechanism [123]. EMRIs can also be made from 
the tidal disruption of binaries that pass close to the MRH [124]. possibly ejecting the 
hypervelority stars observed in our Galaxy (see. e.g., [125]); and from massive-star 
formation and rapid evolution in the MBHs accretion disk [126} . Different mechanisms 
will lead to different EMRi eccentricities and inclinations, evident in the GW signal 
[ 1241 . 

The detection of even a few EMRIs will provide a completely new probe of dense 
stellar systems, characterizing the mechanisms that shape stellar dynamics in the 
galactic nuclei, and recovering information about the MRH, the compact object, and 
the EMRI orbit with unprecedented precision [116]. Especially coveted prizes will 
be accurate masses for LO^IO 7 A/ 0 MBHs in small, non-active galaxies, which will 
shed light on galaxy MRH correlations at the low-mass end; MRH spins, which will 
illuminate the mechanism of MBH growth by mergers and accretion (see Sec. 4); 
as well m stellar- BH masses, which will provide insight on stellar formation in the 
extreme conditions of dense galactic nuclei. The key to measurement precision is 
the fact that the compact object behaves as a test particle in the background MBH 
geometry over hundreds of thousands of relativistic orbits in a year; Hie resulting G W 
radiation encodes the details of both Hie geometry and the orbit [127, 128, 129, 130]. 

To assess the eLISA science performance on EMRIs, we model their very 
complicated signals [131 using the Barack Cutler (BO) phenomenological waveforms 
[132], which are not sufficiently accurate fur detection, but capture the character and 
complexity of EMRI waveforms. We complement this analysis with more realistic 
Tenkohky- based (TB) waveforms obtained by solving the perturbative equations for 
the BH geometry in the presence of the inspiraling body [133]; these have been 
tabulated for circular- equatorial orbits and for some values of MBH spin [130, 134 j. 

To evaluate expected EMRI detection horizons and detection rates, we perform 
a Monte Carlo over 500,000 realizations of the source parameters, taking MRH rest 
mass in [10**, 5 x 10 6 ] AT with a uniform log A/* distribution; MRH spin uniformly 
in [0,0.95]: compact-body mass of 10 A/.,, representative of a stellar-mass RH; orbit 
eccentricity before the fund plunge uniformly in [0.05.0.4]: and all orbital angles and 
phases with the appropriate uniform distributions on the circle or sphere, with an equal 
number of prograde and retrograde orbits. We take the poorly known EMRI formation 
rate to scale with MBH mass ns 400 Gyr l (M */ 3 x l(l <} M . } -019 [135, 130. 137], and 
we distribute systems uniformly in eonioving volume. Our assumptions are consistent 
with the MBH mass function derived from the observed galaxy luminosity function 
using the M* rr relation, and excluding Sc-Sd galaxies [138, 139. 134]. 

We further assume an observation time of two years, consider EMRIs in the last 
live years of their orbit [134], and require a detection HNR = 20 [140, 141. 142], The 
left panel of Fig. Ci shows the resulting nuuimum horizon reds hi ft for BC waveforms, 
as a function of MBH rest mass- that is, it shows the : at which an optimally oriented 
source with the most favorable MBH and orbit parameters (as found in the Monte 
Carlo) achieves the detection SXR, Thus, EMRIs in the eLISA range will be detectable 
jus far z 0.7, By contrast, EM observations of 10 4 10° M© MBHs are possible in 
the local Universe out to : 0.1. The right panel plots the distribution of SNRs as a 
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Figure 6* Loft: maximum horizon redshift vk, MBH n.?st mass. BC waveforms 
{red curve); averaged horizon raiahift va. MBH rest mass, TB waveforms with 
a 9 /hU 0 and 0,9, Assumptions are given in the main text,; the maximum is 
computed as the highest z with SNR > 20 in a given mass bin, Right: maximum 
EMRI SNR vs. redshift, BC waveforms, 

function of z y which shows that nearby EMRIs in the local Universe will yield SNRs 
of many tens. 

For comparison, the left panel of Fig. 6 shows also the horizons computed with 
sky- and orientation-averaged SNRs, using TB waveforms from circular-equatorial 
orbits with MBH spins a, /A/* = 0 ami 0.9. The difference between the BC and TB 
curves is consistent with the effects of sky- averaging: SNRs for optimally oriented 
systems are expected to be 2.5 times higher than averaged SNRs. The a, /A/* 0.9 

systems are favored because high MBH spin allows for orbits closer to the event horizon 
and higher GW frequencies, which shifts the peak eLISA sensitivity to higher masses. 

The resulting number of expected eLISA detections over two years is ^ 50, 
as evaluated with the BC-waveform Monte Carlo, and ^ 30/35/55 (for a 9 /M 9 = 
0/0. 5/0.9), as evaluated with TB- waveform sky-averaged horizons. The higher TB 
event rate is explained by the inclusion of eccentric systems, which radiate more energy 
in the eLISA band, and it should be more reliable because of the broad sampling 
of source parameters. Remember however that EMRI rates are highly uncertain 
[1 16, 135, 136. 123], Even with as few as 10 events, the slope of the MBH mass function 
in the IQ 4 10 6 A/., range can be determined to 0.3, the current level of observational 
uncertainty [143], 

Because EMRI waveforms are such complex and sensitive functions of the 
source parameters, these will be estimated accurately whenever an EMRI is detected 
(140, 141, 142]. In particular, we expect to measure the MBH mass and spin, as well 
as the compact-body mass and eccentricity to better than a part in iO 3 [132], As an 
example. Fig. 7 shows the posterior distributions of the best-determined parameters for 
a z = 0.55 source detected by eLISA with SNR 25, as computed with the Markov 
Chain Monte Carlo algorithm of [144]: for this source, the luminosity distance Di 
would lie determined to 1%, and the sky location to 0.2 deg 2 . Even with relatively low 
SNR, parameter-estimation accuracy is excellent. In general, we find that the eLISA 
and LISA parameter-estimation performance is very similar for EMRIs detected with 
the same SNR (but of course different distances), so the reader can refer tu treatments 
for LISA in the literature 132, 145, 146. 142]. 
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posterior parameter distributims 
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Figure 7. Posterior probability plot for source parameters (MBH rvsi ttt&ss A/*, 
MBH spin a*, com pact- body mass in, and orbit eccentricity at plunge e), in the 
SNR 25 detection of a 10 V 10° A/ 0 EMRI aU = 0.55, with a./A/. = 0.7 and 
Cp lunge : 0.25, 


6* Precision measurements of strong gravity 

Einstein's theory of gravity, general relativity (GR), has been rested rigorously in 
the Solar system and in binary pulsars [117, 14$]; these tests* however, probe only 
tin' weak- field regime where the characteristic perturbative parameter e = v~ jc z ^ 
GM/(Rc 2 ) is very small, ^ 10 -6 li)~ H (here v is the velocity of gravitating bodies, 
M their mass, and /? their separation). By contrast, eLISA’s GW observations 
of coalescing MBHs (Sec* 4} and of EMRIs (Sec* 5} will allow us to confront GR 
with precision measurements of its dynamical, strong-field regime T and to verify that 
astrophysical BHs are really the Kerr mathematical solutions predicted by GR* 

Before considering the GR tests possible with each of these sources, we note 
that* by the second half of this decade, second-generation ground-based detectors are 
expected to routinely observe the coalescences of stellar- mass BHs and (possibly) of 
asymmetric systems such as a NS mspiraling into a 100 M y BH. However, they will do 
so with 10 100 times lower SNRs than eLISA (for the brightest sources)* and for up 
to 1.000 times fewer GW cycles; thus, eLISA will test our understanding of gravity in 
the most extreme conditions with a precision that is two orders of magnitude better 
than that achievable from the ground. (Although most of the references cited in the 
rest of this section were developed for LISA, their broad conclusions are applicable to 
sources detected with comparable SNRs by eLISA.) 

All three phases of MBH coalescence offer opportunities for precision 
measurements. The year-long inspiral signals can be examined for evidence of a 
massive graviton* resulting in a frequency-dependent phase shift of the waveform, 
improving current Solar-system bounds 140 , 150 ]: they can yield stringent constraints 
on ot her theories with deviations from GR parametrized by a set of global parameters, 
such as massless and massive Brans- Dieke theories 151. 152]* theories with an evolving 
gravitational constant [153], Loren tz- violating modifications of GR [154]: last, various 
authors have considered testing inspiral waves for hypothetical, generic modifications 
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of their amplitude and phasing [155, 156, 157, 158]. 

The merger of comparable-mass MBH binaries produces an enormously powerful 
GW burst, which eLISA will measure with SNR as high as n few hundred, even at 
cosmological distances. The MBH masses and spin can be determined with high 
accuracy from the inspiral waveform; given these physical parameters, numerical 
relativity can predict the shape of the merger waveform, as well as the mass and spin 
of the final remnant MBH [159], and these can be compared directly with observations, 
providing an ideal test of pure GR in a highly dynamical, strong-field regime. 

The frequencies and damping times of the quasinortnal modes (QNMs) in the final 
ringdown [160] are completely determined by the mass and the spin of the remnant, 
and therefore can be used to measure them [161, 162], while their relative amplitudes 
hold information about the pre-merger binary [163], again providing a check of 
consistency between GR predictions for the phases of coalescence. Furthermore, 
the measurement oF at least two QNMs [162] will test the Kerr- ness of the MBH 
U64] against exotic proposals such as boson stars and gravastars [165, 166, 167, 168]. 
Modifications of GR that lead to different emission would also he apparent [169, 170]. 

EMRIs are expected to be very dean astrophysical systems, except perhaps 
in few systems with strong interactions with the accretion disk [171, 172, 173], or 
with perturbations due to a second nearby MBH or star [174, 175]* Over day-long 
timescales, EMRI orbits are essentially geodesics of the background geometry; on 
longer timescales, the loss of energy anti angular momentum to GWs causes a slow 
change of the geodesic parameters. In the last few years of their evolution, as observed 
by eLISA, EMRI orbits are highly relativistic (i? < 19/?*) and display extreme forms 
of penastron and orbital plane precession. Indeed, EMRI GWs encode all the mass and 
current mtdtipolesof the MBH [127. 176]* which for a Kerr BH are uniquely determined 
by the mass and spin alone (another manifestation of the “no- hair" theorem). For 
EMRIs with SNR = 30, eLISA will measure mass and spin to a part in 10'*-1()L and 
the mass quadrupole moment M 2 to a part in lQ^-lO 4 t thus testing the no-hair 
theorem directly [129], See [177, 178] for reviews of different ways to test the nature 
of astrophysical BHs. 

Other tests of the Kerr-ness of the central massive object have been proposed: for 
a boson star, the EMRI signal would not shut off after the last stable orbit [179j; for a 
gravastar, QNMs could he excited resonantly [168]; for certain uou-Kerr axisymmetrie 
geometries, orbits could become ergodic or experience resonances [180, 181): for 
“bumpy” BHs, orbits would again carry distinctive signatures [127. 182* 183, 184|, 
Modifications in EMRI GWs would also arise if the true theory of gravity is 
in fact different from GR. as are dynamical Cheni-Situoiis theory [185. 186j« 
scalar tensor theories (with observable effects in NS BH systems where the NS carries 
scalar charge [151. 187)). Randall Sundrum-inspired braneworld models [188, 189], 
theories with axions that give rise to "Boating orbits” [190, 191]. as well as generic, 
phenomenologically parametrized theories [192]. 

7. Cosmology and new physics from the early Universe 

GWs produced after the Big Bang form a fossil radiation: expansion prevents them 
from reaching thermal equilibrium with the other components because of the weakness 
of the gravitational interaction. Thus, relic GWs carry information about the first 
instants of the Universe. If their wavelength is set by the apparent horizon size 
c/H . c(a/d)„ at the time of production, when the temperature of the Universe 
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Figure 8* (Prom [193].) Spectra of stochastic backgrounds from cosmic strings 
for large loops (0.1 the horizon sixe. solid lines), for two values of the string tension 
Gft/c * spanning a range of scenarios motivated by braneworld inflation; and for 
small loops (with size o — 5Q«Gp, dashed tine). The cosmic-string spectrum is 
dislingulshabjy different from that of first-order phase transitions or any other 
predicted source: it has nearly constant energy per logarithmic frequency interval 
over many decades at high frequencies, and Falls off after a peak at low frequencies, 
since large string loops are rare and radiate slowly. Cosmic strings may also 
produce distinctive bursts, produced by a sharply bent bits of string moving at 
nearly the speed of light [212, 213, 21-1, 215]. 


is r„, the redshifted frequency is 

/ « 10- 4 Uz^H~x iSS « 10- 4 Hz . (3) 

so the eLISA frequency band corresponds to the horizon at and beyond the jPetuscale 
frontier* of fundamental physics. This allows eLISA to probe bulk motions at times 
about 3 x 10” 18 3 x 10“ 10 s after the Big Bang* a period not directly accessible 
with any other technique. Taking a typical broad spectrum into account, eLISA has 
the sc ii s i t i v i ty to c ie t.< xl cos n lologic al b a c *kg n ) un d s cat \ sec I ! n r n t* w physi c. s at eue rgie s 
^ 0.1 LOOOTeVj if more than a (modest) fraction ^ 10“ 5 of the energy density is 
converted to GWs at the time of production. 

Various sources of cosmological GW backgrounds are presented in detail in 
[193], They include first-order phase transitions, resulting in bubble micleation and 
growth, and subsequent bubble collisions and turbulence [194, 195, I9G, 107. 108]: 
the dynamics of stabilisation for the extra dimensions required by superstring theory 
199, 200]. which may also appear as non-Newtonian gravity in laboratory experiments 
at the sub-nun scale: networks of cosmic (super-) strings [201. 202], which continuously 
produce loops that decay into GWs {see Fig. 8): the transition between inflation 
and the hot Big Bang in the process of preheating [203, 204, 205, 206, 207 : and 
the amplification of quantum vacuum fluctuations in some unconventional versions of 
inflation (208. 209. 210 . Although the two-arm eLISA does not provide a Sagnac TDI 
combination 211] to calibrate instrument noise against possible GW backgrounds, 
the clear spectral dependence predicted for some of these phenomena provides an 
observational handle* as long as the background lies above the eLISA sensitivity curve. 
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As discussed in Set;, 4, observations of GWs from MBH binaries probe the 
assembly of cosmic structures. In addition, binaries can serve as standard sirens 
to measure cosmological parameters [216, 217] because, as discussed around Eq. 
(2), measuring the amplitude and frequency evolution of a binary signal yields the 
absolute luminosity distance to the source. However, binary GWs cannot provide 
the source’s red shift unless the other source parameters are known independently 
(because the rest mass of the binary is the only length/ time scale in the waveform, 
the frequency evolution of a redshifted signal is indistinguishable from the signal from 
a heavier binary). The optical redshift of the host galaxy can be obtained if an EM 
counterpart to MBH coalescence is observed (see, e.g., (218, 219, 220], and {221] for a 
recent review}. 

While there are many uncertainties in the nature and strength of such 
counterparts, some may be observable in the local Universe, At z < 1, we expect that 
eLISA MBH-inspiral measurements could provide sky locations to better than 400 deg 2 
for 50% of sources, and to 10 deg 2 for 11%, (The inclusion of merger and ringdown in 
the analysis should further improve these numbers.) Such large areas will be covered 
frequently and deeply by optical and radio surveys such as LSST [222] and the VAST 
project [223], identifying sufficiently distinctive transients. The accurate knowledge 
of the counterpart’s redshift and position would then improve the uncertainty of 
GW- deter mined parameters, with known to 1% for 60% of sources, and 5% for 
87%, Such precise luminosity distance-redshift measurements will be complementary 
to other cosmographical campaigns [224, 225], and wall improve the estimation of 
cosmological parameters. Even without counterparts, one may proceed by considering 
all possible hosts in a distance position error box, and enforcing consistency between 
multiple GW events [226]: this should be possible for MBH binaries (and EMRIs [227]} 
in the local Universe, yielding the Hubble constant to a few percent. 


8. Conclusions 

While LISA was always meant to be the definitive mission in its frequency band, eLISA 
is being designed to provide the maximum science within a cost cap. Nevertheless, 
as described above, eLISA will achieve a great, part of the LISA science goals. It will 
represent the culmination of twenty years of exciting, painstaking work, pioneering 
the new science of observational low-frequency GW astronomy. It will t ruly begin to 
unveil the hidden, distant Universe. May it Hy soon, and safe. 
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